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Abstract 

We consider the integral of fractional Brownian motion (IFBM) and its 
functionals £t on the intervals (0, T) and (— T, T) of the following types: 
the maximum Mt, the position of the maximum, the occupation time 



oo, 



<**] 1 above zero etc. We show how the asymptotics of P(£t < 1) = pr, T 

is related to the Hausdorff dimension of Lagrangian regular points for 
the inviscid Burgers equation with FBM initial velocity. We produce 
computational evidence in favor of a power asymptotics for pr- The data 
do not reject the hypothesis that the exponent 6 of the power law is related 
to the similarity parameter H of fractional Brownian motion as follows: 
6 = -(1 - H) for the interval (-T,T) and 9 = -H(l - H) for (0,T). The 
point is special in that IFBM and its derivative both vanish there. 
OO ■ Key words: Fractional Brownian motion, Burgers equation, fractal- 

' ity, large excursions 

(N 

O 
^— > 

s 

> 

X 
S3 



*Observatoire de la Cote d'Azur, CNRS UMR 6529, BPU 229, 06304, Nice Cedex 4, France, 
molchan@mitp.ru 

t International Institute of Earthquake Prediction Theory and Mathematical Geophysics 79, 
b2, Warshavskoe shosse 113556 Moscow, Russia, khokhlov@mitp.ru, khokhlov@ipgp.jussieu.fr 



1 



1. Introduction 

Sinai [18] and Frisch and associates [17] initiated in 1992 the study of frac- 
tal and multifractal properties of solutions of the inviscid Burgers equation 
with initial velocity uq(x) specified by a self-similar random process. That last 
circumstance guarantees that the solution is self-similar in the large. In par- 
ticular, one could be interested in finding the Hausdorff dimension of the set 
of regular Lagrangian points S that describe the initial locations of those fluid 
particles which have not collided until a fixed time to- The original model of 
uq(x) was fractional Brownian motion (FBM), bn(x), with similarity parameter 
< H < 1. 

By now the Sinai- Frisch program has been carried out for special Markovian 
models of uq(x) alone: Sinai [18] has found the dimension S for Brownian motion 
case, i.e., uo(x) — bi/ 2 (x); Bertoin [3] discovered for this case that the solution 
u{t = to,x) admits of an exact probabilistic description in terms of a stable Levy 
process. One can then find a multifractal description of the solution x — > u(to, x) 
(the relevant references are [6, 7, 8]). Additionally, Bertoin [3] found that the 
Hausdorff dimension of Lagrangian regular points is h, if Uo(x) is a stable Levy 
process of index a = h^ 1 e (1, 2] with no positive jumps (see also [20]). 

The nonmarkovian case Uq(x) = bjj(x), H ^ 1/2 has proved extremely diffi- 
cult for analysis. Handa [5] found simple arguments to derive a lower bound on 
the dimension of 5, namely, dimS* > H. The exact equality dim5 = H is known 
as a hypothesis [17, 19] since 1992. Among methods developed for analyzing the 
nonmarkovian case uo(x), the Sinai approach is of particular interest. For the 
case uo(x) = 61/2(2;) this method [18] relates the estimation of the dimension of 
S to the asymptotic behavior of the probability 

PT = P{£(x) < l, « A T } (1) 
for integral Brownian motion 

£(a;)= / b 1/2 (s)ds and A T = (0,T), T > 1. (2) 
Jo 

As a matter of fact (see below) , one has to deal with a problem that is rather 
popular in physical and technical applications: find the probability of a large 
excursion for a random process rj(x), i.e., P{rj(x) > 0, 1 < x < T}, T> 1. A 
review of the problem can be found in [14]. Sinai has shown that the quantity 
Pt ■ T 1 / 4 is bounded away from and 00 as T — > 00 under the conditions (2). 
That estimate was repeatedly refined and generalized [9, 10, 11]. 

We show below that the upper bound dim5 < H under the conditions 
uq(x) = bu{x) follows from an estimate of pt for the integral of fractional 
Brownian motion (IFBM): £(x) = J Q X bn(s)ds when considered in the bilaterally 
expanding interval At = (—T,T). 
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The work [15] clarifies the asymptotic problem of px for intervals (0, T) and 
(— T, T) in the case of fractional Brownian motion: £(x) = bn{x). It transpires 
that in this case 



lnp T = -(l-i2")mT(l + o(l)), A T = 



(0,T). 



On the other hand, when At = {x : \x\ < T}, the leading term in the log 
asymptotics of pt is independent of H. More generally, suppose bn(x), x € R d 
is FBM with multidimensional time; in that case 



The last asymptotics is due to the fact that the probability density for the 
position of the maximum of FBM exists in the sphere {|x| < 1}. A generalization 
of this fact is given below. 

We present theoretical and computational evidence in favor of the following 
asymptotics for IFBM: 



The first of these asymptotic expressions corroborates the hypothesis dimS 1 = 
iJ, so is not unexpected, while the second is, considering that the exponent 
0{H) = H(l - H) has the point of symmetry H = 1/2. 

Because IFBM is a self-similar process, the distribution of its maximum in 
A = (0,1) or (-1,1), F max (x), is related to p T through p T = F max (T-< 1+ ^). 
Importantly, our calculation was performed for a series of statistics: the maxi- 
mum M = max IFBM; the position of the maximum M, \G\; the occupation 

time A + = j A l^^yodx of IFBM above zero; and the rightmost zero of IFBM 
in (0,T), Z. The distributions of these statistics (one should use F max (x 1+H ) 
when M is considered) have identical asymptotics as x — > 0, but depend on in- 
terval type: A = (0, 1) or (-1,1). When A = (—1, 1), they provide independent 
evidence in favor of the hypothesis dimS* = H. 

The rest of this paper is organized as follows. Section 2 reduces the evalua- 
tion of dimS 1 < H to the asymptotic distributions of M, G, A + and Z near zero. 
Section 3 discusses the modeling of IFBM, while section 4 presents numerical 
evaluations of the distributions listed above and some theoretical arguments to 
support our conclusions. 

2. Regular Lagrangian points and the nonexceedance of level 

We now define more exactly the notions used in Introduction. We consider 
the Burgers equation 



In Pt 



dlnT(l + o(l)), A T = {x e R : \x\ < T}. 



In pt = 



-(l-*01nT(l + o(l)), A r = (-T,T) 
-fl-(l-ff)InT(l + o(l)), A r = (0,T). 
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d t u + ud x u = vu xx , v i 0, 



(3) 



with continuous initial conditions u(0,x) = u (x) and the velocity potential 
U(x) = L uo(x)dx = o(x 2 ), x — » oo. The solution at to = 1 has the form 
u(x) = x — a(x), where a(x) can be found from U{x) as follows. Construct 
a convex minorant C(x) for U{x) + x 2 /2. In that case its derivative C'{x) 
is nondecreasing and has finite limits from the left and from the right. We 
now complete the definition of C'{x) in continuity on the right. In that case, 
according to Hopf (see, e.g., [18, 21]), a(x) is identical with the inverse func- 
tion of C'(x). The set of points where C'(x) is increasing, i.e., the topological 
support of the measure dC'(x) or the closure of the set {a(x),x € R}, de- 
fines the set S of regular Langrangian points in the Burgers problem. The 
dynamics of completely inelastic particles on R 1 can be related to the Burgers 
equation: each infinitesimal particle located at x has a mass dx and an initial 
momentum dll(x). On colliding the particles coalesce and continue movement 
following the conservation laws of mass and momentum. Particles that have 
not collided until time to = 1 make up the set S in the Lagrangian coordi- 
nates. The initial conditions uq(x) will be considered to be fractional Brownian 
motion bu{x) 1 i.e., a Gaussian process with zero mean and structural function 
E\bn(x) — bH{y)\ 2 = \x — y\ 2H where < H < 1. In virtue of the Kolmogorov 
theorem the paths of bn{x) can be treated as continuous a.s. The process 

bn{x) is self-similar, i.e., 6#(Ax) = A ff 6#(x), where = denotes equality of 
finite-dimensional distributions. 

Theorem 1. 1. The set of regular Lagrangian points S in the Burgers 
problem (3) with uq{x) = bn(x) has a.s. dimension H, if for any e > and 
T — > oo one of the following requirements is fulfilled: 

(A) P{y(x) := f b H (s)ds < 1, x e A T } < T^ 1 ""^, 

Jo 

(B) P(y(x) < 0, x e A T , \x\ > 1) < T^ 1 -^, 
(CO P(\G(A T )\ < 1) <T - (1 - H '> +e 7 

(D) P{( l y{x)>0 dx<l, |G(A T )| <T}<T-( 1 - ff )+ £ , 

where At = (— T, T), G(At) is the position of the maximum of y(x) in At- 

2. If one of type A — D probabilities pr has an asymptotics of the form 
\ogPT = —0\ogT(l + o(l)), the probabilities of the other types have the same 
asymptotics. This statement also holds for A T = (0, T) with the probability 
P(Z T < 1) in addition to (A — D), where Z T is the rightmost zero of y(x) in 
(0,T). 

The proof of the theorem will be preceded by two lemmas. 
Lemma 1. dim5 < H, if for any e > there exists a So — So(e) such that 
one has for arbitrary x e R 1 : 
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P(S n B(x, 6)^<j))< S^- H ^ E , 6<6 a 



(4) 



where B(x, S) is a ball of radius 6 centered at x. 

Proof. Cover the interval A = [a, b] with intervals Bi{8) of length 6 with 
overlappings of length 5/2. Consider the measure n(dx) = dC'(x) with support 
S. The elements Bi in {Bi(S)} for which fi(Bi) > will then form a cover Sn A. 
In view of (4) 

Ej2m8)\ H+2s = E^lB^+^l^^ 
< 5 H+2e ■ 2\^\5- 1 ■ St 1 - 11 ^ = cd E 

where |A| is the length of A. By Chebyshev's inequality 

P(^\Bi\ H+2e >a)< c6 E /a. 

Consider a sequence 8 n such that ^ S„ < 00 ■ The Borel-Cantelli lemma then 
yields 

]T \Bi(5 n )\ H+2s < a, n>n(oj) 
Since a is arbitrary: 

limsup^ \Bi{5 n )\ H+2e = a.s. 

n 

However, in that case one has dim(5 n A) < H + 2e. Since e > and A are 
arbitrary one has dimS* < H. 

Lemma 2. The conditions of Lemma 1 are fulfilled, if 

P{ [ b H (s)ds < 1, |z| <T}< T-( 1 - Ef )+ £ , Ve > 0. 
Jo 

as T — > oo. 

Proof. The process y(x) = f£ bH(s)ds + x 2 /2 can be represented in the form 

y(x) = [ (b H (s) + s)ds + (b H (c) + c)(x-c) + 
Jo 

[b H {s) - b H {c) + (s - c)]ds = L(x') + J (b H (s) + s)ds, 
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where L(x') is a linear function of x' = x — c, and b H (x) = b H (c + x) — b H (c) = 

bn(x). The convex minorants of y and y = J Q X (6ij(s) + s)ds differ by the linear 
function L(x'). Hence the fractal properties of the measure fi(dx) = dC'(x) 
are invariant under translation along the x-axis (this observation is due to U. 
Frisch). Consequently, it is sufficient to prove (4) for S' = S fl (—5/2,5/2). 

Let A = (—<5/2, 5/2) contain a point of growth xq for the measure dfi. That 
means that the curve f(x) — U(x) + x 2 /2 and its convex minorant C(x) do not 
lie below the tangent of f(x) at the point xo, and C(xq) = f(xo). The event 
{x € A}, to be called A here, can be written as 

A = {3x Q :\x \ < 5/2; [ (b H (s) + s)ds 

Jo 

> / (b H (s) + s)ds + (b H (x ) + x Q )(x - x Q ), Vx £ R 1 }. 
Jo 

Let us modify event A to become A\, i.e., we assume that the equality in 
the formulation of A is true for < 1 only. To emphasize the fact that A\ 
depends on the process bn(x) + x = £(x), we will write A\ = A\ [£]. 

One has 



P(A) < P{A,) - E l AllbH+v] = E l Al[5a] 7r(M, (5) 

where ip(x) — x and ir is the Radon-Nikodim derivative of two Gaussian mea- 
sures corresponding to the processes bn — f and bn in [-1,1]. Note that bn is an 
FBM process. The function ip is smooth and vanishes at zero. For this reason 
the above measures are mutually absolutely continuous [16]. By the Cameron- 
Martin relation ln7r(&ir) is a Gaussian variable with mean —c? H /2 and variance 
c H , where ch = |M| and || • || is the norm in Hilbert space Hb of functions on 
A = [—1, 1] with reproducing kernel B(x,y) = EbH{x)bu{y)- The constant ch 
is finite and can be found in explicit form as indicated by Molchan and Golosov 
[16]. 

Applying Holder's inequality to the right-hand side of (5), one gets 

P(A) < P(A 1 [b H ]) 1 - e (Ew 1 ' e ) e = P(A 1 [b H ]) 1 -"c £ , (6) 

where c £ — exp(i(e _1 — l)c 2 H ). 

We now evaluate P(Ai[bn])- One has 

P(A![b H ]) = P{3x :\x \<5/2; 

fX rXo 

/ b H (s)ds > / b H {s)ds + b H {xo)(x - x ), \x\ < 1} 
Jo Jo 

= P{3x :\x \<l/2;[ b H {s)ds > a{x ) + b(x )x, \x\ < T}, 
Jo 
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where 



rxo 

T = 5-\ \a(x )\ = 1/ b H (s)ds-b H (x Q )x \ < 2M, 

Jo 

\b(x )\ = \b H (x )\ <M= max \b H (x)\. 
We will use the Fernique inequality [4] 

P(M > c H u) < exp(-u 2 /2) = T~ a , u > u , 



where u = Ut = \/2a InT, ch being a constant; the value of a will be chosen 
later on. From this it follows that 



P(Ai[6 ff ]) < P{Ai[6 ff ] ) M<c H u T } + T- a 

< P{ / fcff (s)ds > -2c ff u T -c ff u T |x|, |x| < T} + T- a 
Jo 

= P{[ b H {s)ds < u T c H (2 + \x\), \x\ < T} + T- a 
Jo 

= P{[ b H (s)ds<4X^ 1 + 2\x\, \x\ <T'}+T~ a , (7) 
Jo 

where V = T/X T , u T c H = 2A^, At = const • {InT) 1 / 211 . Here we have used 
the fact that bn(x) is a self-similar process and modified the interval |x| < T to 
become |x| < T'. 

Define the function 

= 2x 1 M<1 + 2sgn(x)l| x | >1 = A y [ e «* _ i]!EL* dA . (8) 
In that case (7) can be continued to get 

P{M[b H \) <P{[ (b H (s) - Vl {8))ds < F(x), \x\ < T'} + T~ a , (9) 
Jo 



where 



F(x) 



-x 2 + 2\x\ + 4A T "\ |x| < 1, 

1+4:X T \ |x| > 1. 



When T is large, one has F(x) < 2. For this reason the last estimate will merely 
become less precise, when F is replaced with F(x) = 2. The right-hand side 
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of (9) can be evaluated by repeating the steps that have led to (5, 6). The 
substitution of bu for bn — fi combined with Holder's inequality yield 

PiA^bn}) <P{[ b H (s)ds < 2, \x\ < ry- e cl + T- a , 
Jo 

where cj = cxp((e _1 — 1)6^/2), b 2 H = \\fi\\^ < 1 1 V^i 1 1 So - Here || • \\t is the norm 
on Hb for the interval (— T, T). The spectral representations of the kernel 

B(t, s) = Eb H (t)b H (s) = [ (e" A - l)(e" te A - l)^!" 1 - 2 * dX 
and ipi (see (8)) yield 

yi\\l = kH J 

where k H = J \e tX - 1 1 2 ! A| ~ x ~ 2ir rfA. 
The final result is 

P(A) < PiA^bH}) 1 -^, < ip 1 T 7 e c 1 £ +T- a ) 1 - £ c e , (10) 

where p T ' = P{f* b H (s)ds < 2, \x\ < T'}. 

Let pt < T~^~ H > +£1 for large T. Take a > 1 — H and choose e from the 
requirement c\ ■ c £ — T' ei , i.e., e — cs^ 1 / 'InT' '. Inequality (10) can then be 
continued: 

P(A) < Cl T'- (1 - H ^\ 

where a = exp(2(l - i/)ce^ 1 ). Recalling that T = T(lnT)-^ 2H ■ c 2 , one 
obtains the desired estimate P(A) < T-( 1 - ff )+ 3£l , T > T (si), T = 5' 1 . 

Proof of Theorem 1. The inequality dimS > H was derived by Handa [5]. 
The opposite inequality dimS 1 < H follows from Lemmas 1 and 2 and condition 
(A) of Theorem 1. To prove the theorem under condition (B), we note that the 
event {Jq bH(s)ds < c, \x\ < T} can be represented as 

{ f (b H (s)-cp(s))ds<iP(x), \x\<T}, 
Jo 

where ip, ip are smooth finite functions: ip = when |x| > 1 and ip > when 
| a; | < 1, while ^> = when |x| < 1/2 and \x\ > 1. Repeating the translation 
procedure for the samples: bn(s) — <p(s) — ► 1>h(s) and using Holder's inequality 
we get 



2 sin A 



ttA 2 



|A| 



1 + 2H 



dX < oo, 
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PT ■ = P{[ b H (s)ds < 1, |z| < T} 
Jo 

< c £ P{ [ b H (s)ds < \x\ < T} 1 - 5 

Jo 

< c £ P{ [ b H (s)ds < 0, 1< \x\ < T} 1 - 5 , 

Jo 

where c e = exp(ie _1 • c 2 ), c v < kjf J \ip(\)\ 2 \A\ 1+2H and (p is the Fourier 
transform of ip. One has c v < oo, because y> is smooth and finite. Choose 
e = £t from the requirement c e = £t, where Ct is a slowly varying function. 
Take Ct = hiT, say then e^ 1 = clnlnT. The result is 

f x 

Pt<C t P{ b H (s)ds < 0, 1< \x\ < T)^ £T . 
Jo 

When (B) holds, one has 



Pt 



< (T~( 1 ~ H ^ +£l ) 1 ~ ST £ T < T- {1 - H)+e2 , T> 1, 



i.e., the implication (B) — > (A) is true. 
The inequality 

P( / < 0, 1< |ar| < T) < P(|G(-T,T)| < 1), 

Jo 

where G(A) is the position of the maximum of IFBM in A, yields the implication 
(C) — > (B). Lastly, under (B) the position of the maximum of IFBM is \G\ < 1, 
while the occupation time of IFBM above is below 2. Hence (D) — > (£>). 

Let us prove (A) — ► (C). Below, Gt is the position of the maximum of 
IFBM in (— T, T) and M a = max IFBM. One has 

|a;|<a 

P(\G T \ < 1) < P(\G T \ < 1), Mi < c T ) + P{M l > c T ). 



If ct = V2a lnT, then the Fernique estimate [4] yields 

P{M 1 > c T ) < cT- a ', T > T , 
where a' = a/a 2 and a 2 = £|IFBM(1)| 2 = (2H + 2)" 1 . Also, 

P(|G T | < 1, Mi < c T ) < P(M T < c T ). 
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Since IFBM is self-similar, one has M T = \ 1+H M T ,, when T = XT'. Take 
A from the requirement X 1+H = c T . Then P(M T < c T ) = P(M T > < 1). 
To sum up, 

P(\G T \ < 1) < P(M T , < 1) + o(T- a ), 

where T 1 = cT(lnT)~ p , p = (1 — H)^ 1 , while the parameter a > is arbitrary. 
When a > (1 — H), the implication (A)— >(C) is obvious. 

We are going to prove (A)— >(D). Let be the occupation time of y(x) = 
IFBM(x) above zero in A T = (— T, T). One has 

P{A+ < 1, \G T \ <T)< P(M t < or) + P{M T > c T , A+ < 1, \G T \ < T), 

where ct will be specified below. 

Let At = UA&, A^ = (k, k + 1) and Mfc = max{y(x), x e A^}. If the event 
B = {Mt > ct, A^, < 1, \Gt\ < T} occurs, one will have the following for the 
interval A& which contains Gt- Mk > Ct, y(x) and y'(x) = bu{x) have zeroes 
in Afc. Indeed, if y{x) ^ 0, then y(x) > in A fe and A^ > 1. Consequently, 

P(B) < P{max((y(a;i) - y(x 2 )) > c T , x x ,x 2 € A fe ), S^} := J^Pfc, 
fe fe 

where Sfe means that bn(x) has a zero in A&. 
We are going to evaluate pk- 

p k < P{max[(y(x 1 ) - y(x 2 )), x x ,x 2 € A fe ] > c T , \bH(k)\ < c T /2} 
+ P{\b H (k)\ > or/2, S k } := p kA + Pk , 2 . 

One has 

Pk.2 < P{m&x(\b H {xi) ~ b H (x 2 )\, x r ,x 2 e A fe ) > c T /2} 
= P{max(|fe ff (a; 1 ) - 6^(^2)1, Zi,^ & A ) > c T /2}. 

Here we have used the fact that has stationary increments. In virtue of 

the Fernique inequality [4] 

Pk. 2 < c exp(-i(c T /c 6 ) 2 ), 

where c is an absolute constant, while c& is a function of H. 
One proceeds in a similar manner to evaluate pk.Y- 
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y(xi) - y(x 2 ) = / (6i?(s) - b H (k))dx + b H {k){x 1 - x 2 ). 

JX2 

If max.[y(x\) — y(x 2 )] > ct in A& x A& and |&Ef (fc)l < ct/2, then 

/■Xl 

max / [&ff(s) — &ir(fc)]ds > Ct/2. 

Consequently, 

/■Xl 

Pk,i < P{max[ / (b H (s) - b H (k))ds, x\,x 2 € A fe ] > c T /2} 

/■X2 

= P{max[ / bH{s)ds, xi,x 2 € A ] > ct/2}. 

■/ Xl 

Here again, we have used the relation b H (x) — b H (k) = b H (x — k) with a fixed 
k. The use of the Fernique inequality [4] yields 

Pk,i < c exp(-i(c T /2c s/ ) 2 ) 

where c y is a function of H. Combining the estimates of pk,i and pk, 2 and 
assuming ct = max(o,, c lc y )\J2a logT, one gets 

Pk =Pk,i + Pk,2 < cT~ a . 
However, one then has P{B) < 2cT~ a+1 and 

P(A+ < 1, \G T \ < T) < P(M T < c T ) + 0(T- a+1 ) 
= P{M T , < 1) + 0(T- a+1 ), 

where T = cT(logT)-P, p=(l- if)" 1 . Hence (A) -» (£>). 

Consider the second part of the theorem. Let be the probabilities that 

appear in the first part of the theorem, where denotes the condition A, B, C 
or D. It has been shown above that, when T ^> 1, 



Pt(A) < pr^Bf-^Cr, pt(B)<pt(C), 

p T (C) < PT ,(A)+0(T- a ), PT (D)<p T ,(A) + 0(T- a+1 ), 

where a > is any fixed number, Ct is a slowly varying function of T, St = o(l), 
T -> oo and T = cT{\nT)~P, p = (1 - H)' 1 . A trivial corollary of these is that 
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all the Pt(Q) have the asymptotics lnpx(@) = — 01nT(l + o(l)), provided the 
asymptotics is true for at least a single quantity of the type 6 = A, C or D. 

Our proof has not relied significantly on the type of the interval At- {—T, T) 
or (0,T). For this reason our conclusion that the asymptotics of pr(O) are 
identical also holds for (0, T). 

We conclude by noting that, if A T = (0,T), then p T {B) = P{y{x) > 0, 1 < 
t <T). Consequently if Zt is the rightmost zero of y(x) in (0, T), then P(Zt < 
1) = 2 PT (B). 

3. The generation of IFBM 

We are going to use Monte Carlo techniques in order to evaluate the proba- 
bilities pt(@) with = A, C, D in Theorem 1 for the process y(x) = J Q bH(s)ds 
in the following intervals of At ■ (0,T) and (— T, T). The probabilities in ques- 
tion are small, pr — ► as T f oo, hence the IFBM generation should be exact for 
a discrete sequence {xk,k — 1, ...,T}. Since y(x) is a self-similar process, it is 

sufficient to use integer points Xk = k. In that case {y(k/T)} = {T~( 1+H ^y{k)}, 
while the probabilities can obviously be expressed in terms of the statis- 

tics M — maxj/(x), G = argmaxy(x), A + = J A l y>a dx of the process 

Ai Ai 1 

{y{x), x e Ai}, where Ai = (0, 1) or (—1, 1), as follows: 

p T {A) - F M (T-^- H ^); PT (C) = F^T- 1 ); p T (D) - ^^"^^(l - 0) 

where F% is the distribution of ^ and Fa is the conditional distribution of A + 
given \G\ ^ 1. 

The generation of {y(k), k = 0, T}. The sequence {y(k), k = 0, T} is 
Gaussian and has stationary second increments, i.e., the sequence 

Vk=y(k-l)-2y(k) + y(k + l), fc = l,...,T-l, (11) 
has a Toeplitz correlation matrix [l-ii-j], where 

Aife = c q [\k - 2\" - 4|fc - l\ q + 6|fc|« - 4\k + 1\" + \k + 2\«] (12) 

and c q = [2q(q - l)]" 1 , q = 2H + 2. 

The second differences (11) combined with the initial conditions y(0) = 
and y(l) are sufficient to uniquely reconstruct the sequence {y(k),k = 0, T}. 
One can assign y(l) by using the decomposion y(l) = y(l) + y (1) into the 
predictable part y(l) = E(y(l)\r]i, rfr-i) of y(l) and the part y _L (l) that 
cannot be predicted from the data {r/i, i — 1, T — 1)}. In that case 

T-l 

y(l) = ^2 ZkVk + o-eo (13) 
fe=i 
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Hero, z = (z\, Zt-i)' is the solution of the linear equation: 

\Mi-AT~ 1 * = m > ( 14 ) 
where the vector m has the components 

m k = Ey(l) Vk = Afelfcl^ 1 - \k\" + \k - 1|% 

and A is the difference operator of second order: A/(fc) = f(k — 1) — 2f(k) + 
f(k + 1). The second term is <reo = y _L (l), where e is the standard Gaussian 
random variable which is independent of {r)\, 7?t-i}; 

a 2 = E[y^{\)] 2 = q- 1 - £ z k m kl 

l<k<T 

because q^ 1 = E\y(l)\ 2 and E\y(l)\ 2 = Y^ 1 z km k . 

It thus appears that the exact generation of the sequence {y(k), k — 0, T} 
reduces to the generation of the stationary Gaussian sequence {rj kl k = 1, ...,T— 
1)} with correlation function (12) and to the solution of the linear equation (14). 

The generation of {y{k), \k\ < T/2}. For generating y(x) in a bilateral 
interval, we note the following. Assume that y(x) is IFBM in (0,T), while y(x) 
is IFBM in (-T', T"),T' + T" = T. In that case 

{y(x) y{T>) y'(T')(x — T'), x e (0, T)} 4 - T'), x e (0, T)} 

The left-hand side provides a key to how one is to transform the sequence 
{y(k), k = 0, ...,T} into an IFBM sequence that starts from the point < k < 
T. To do this one must also find the derivative y'(ko). In a similar way as above: 

y'(k ) = E{y'(k )\vu -,Vt-i} + E{y'(k a )\e Q } + y'^{k ) 

where the first two terms correspond to the predictable part of y'(ko) based on 
the data {rji, ...,t]t-i,So}, while the third term corresponds to the unpredictable 
part of y'(ko). The predictable part is 

T-l 

£{y(MI»?i, ...,77t-i,£o} = z 'kVk + ae , 

k=i 

where (z[, z^-i) is t ne solution of (14) with the right-hand side m' = (m' 1; ...,m' T _ 1 ). 
The components of m' are 

m' k = Ey'{k )n k = + \k - fc| <? - 1 sgn(fc - k))qc q . 
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One has a = Ey'(ko)so- From (13) one derives Ey(l)y' (k ) = o-a+Y^k=i z kfn'k- 
Hence 

T-1 

aa = c q - q[(q - l)fc<f 2 + 1 - fc^ 1 + (k - l)^ 1 } - ]T z k m' k . 

fc=i 

One has y' ± (ko) = cr'e' where e' is a standard Gaussian variable that is 
independent of (rfi, r? n -i, e). The variance of the unpredictable part cr' 2 can 
be found from the relation 

fc ^^£;[ y '(fc )] 2 -E« + « 2 + ^ 2 - 
fe=i 

To sum up, the exact generation of {y(k), k = 0, T; y(fc ) = y'{k ) = 0} 
requires that an equation like (14) should be solved twice. 

The generation of {rjk}- Bardet et al. [2] provide a review of the meth- 
ods which allow generation of Gaussian stationary sequences with a prescribed 
correlation function. We use the progressive Schur algorithm [1], which is a 
Levinson-Durbin method. The Generalized Schur algorithm can be used in the 
framework of this method for fast solution of equations like (14) by the Gohberg- 
Semenkul formula [1]. The generation of {r]k, k = 1, T — 1} by this method 
requires 0(T 2 ) operations. The computation is organized so as to minimize the 
amount of calculation needed for generating N IFBM samples; the computa- 
tional complexity is a linear function of N and the storage capacity is of order 
0(T 2 ). 

The parameters T and N are equally important in the problem discussed. 
However, it appears that the increase of T would not be effective from a pri- 
ori considerations. The argument is as follows. We are interested in the dis- 
tributions of M, G, A+ near zero where they are expected to behave like 
x e £(x), where £ is a slowly varying function. Since the approximation to 
IFBM is discrete, these distributions contain a positive atom at zero of size 
p (T) = P(y(k) < 0, k e A T ). The probability is doubled for the statistic Z 
(the rightmost zero of y(t),t € (0, T)). 

Note that po(T) > P(y(x), \x\ > 1, x <G At). Therefore, if the bounds given 
by Theorem 1 are explicit, Po{T) should not decrease faster than T~ e £ (T). The 
rate of decrease becomes very low, when 9 < 1. The expected value is 9 = 1 — H 
for the interval (— T, T) and 9 < 1/4 for (0,T) (see below). 

Given the above situation, it only remains to increase N alone. This will 
provide good accuracy for empirical estimates of the distributions of M, G, A + , 
Z for discrete time. The discrete-time distributions cannot converge to the basic 
continuous ones faster that T~ e ( H ) that we have at x = 0. 

4. Evaluation of 9(H) and related results 
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The bilateral IFBM process. 

Figure 1 presents estimates of the distribution of \G\/ 2 \ (G a is the position of 
the maximum of y(x) = 6_y(s)ds in the interval (—a, a)). The distributions 
are given for H : 0.1 — 0.9 at increments of 0.1. The estimates are based on 
N = 50,000 samples of y(x) with sample size T = 8194. We use \Gi/ 2 \ to 
demonstrate the asymptotics of pr from Theorem 1, because in this case the 
distribution discontinuity at x — bends the graph of the distribution on a 
log-log scale near zero to a lesser extent. We recall that the discontinuity of size 
Po(T) occurs for all the statistics M, G and A + , and is due to the presence of 
negative excursions in y(x) on a discrete grid of x. The curves in Fig. 1 are well 
consistent with the asymptotics of type Fq(x) = x e £(x), where £ is a slowly 
varying function and 9 = 1 — H. Assuming C(x) = constant in (x-,x + ), one 
can construct the estimate of maximum likelihood 9 for 9. 
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Figure 1: Distributions of \G\ for the position of maximum G of the IFBM 
process in interval (— |, |). Various values of H are shown at the left hand part 
of each curve. 

The estimate 9 was computed in intervals of the form (10~ 3 , 10~ 2 ) x i, i = 
1 — 5. The choice of the initial point 0.001 is related to the fact that the sampling 
interval over time is A = T _1 = 0.00012, so that all curves in (A, 10A) slightly 
change their slopes due to the discontinuity of discrete-time distributions at zero. 
The deviation of 9 from the hypothetical 6*o = 1 — H does not exceed 0.03, i.e., 
is less than 6% when H < 0.5. The error is large (> 10%) for 0.6 < H < 0.9 
because of smallness of 9$. For H > 0.6 we have a slow convergence of the 
discrete-time distributions with A = T^ 1 J, and difficulties in the choice of 
small x to estimate 9. 

The argument to be given below provides a partial explanation of the nature 
of the asymptotics 

P(\G T \ < 1) =T-( 1 - ff )£ T . (15) 
The function y{x) is differentiable, hence the position of the global maximum, 
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Gt, belongs to the zero set of buix) or to the end-points of (— T, T). For this 
reason it should seem that the local time l(x) — lim J l\b H ( s )\<eds is the 

natural time scale in our problem of the maximum of y{x), i.e., it is more natural 
to study i/(7) = y(x{l)) instead of y(x), where x(l) is the inverse function of l(x) 
which is continuous on the right. The process was first treated by Vergassola 
et al. [19] and independently used by Isozaki and Watanabe [11] to prove the 
Sinai asymptotics for H = 1/2. It is a known fact [12] that l(x) is a continuous 
self-similar process with parameter h = 1 — H. Consequently, l(T) — 0(T 1 ^ H ), 
and (15) means, roughly speaking, that P(|Gl| < 1) = where Gl is the 

location of the maximum of y(l) in (—L,L). A more exact statement can be 
made. Let Cr,i be slowly varying functions that decrease as T — ► oo. 
Statement 2. If a) Gt and l(±T) are weakly dependent, i.e. 

P{|G T | < 1, |Z(±T)| > T 1 -"^)} > P{\G T \ < 1)£t,2 (16) 

and b) P(\G L \ < 1) < [LC L ]~\ then P(|G T | < 1) < [T^-^Ct.s]' 1 - 
Proof. By (16) we have 

P(|G T | < 1) < C T ] 2 P{\G T \ < 1, \l(±T)\ > T*- h Ct, u 

\l(±G T )\ < (a\nT) H } + C T ^ 2 P{\l{±l)\ > (a\nT) H } (17) 

The second term on the right-hand side admits of an upper bound according to 
[22]: C T \ ■ T- ca where c is an absolute constant. It follows that (17) can be 
continued: 

< L T \ [P{\G L(T) \ < (alnT) H } + 0(T~ ac )], 

where L(T) = T 1 " iJ L T ,i- 

The process y(l) is self-similar; therefore, 

P{\G L[T) \ < (a\nT) H )} = P{\G L .\ < 1} < [L*C L ,}~\ 

where L* — L(T)(a \nT)~ H . Combining the resulting estimates and bearing in 
mind that the constant "a" is arbitrary, one gets Statement 2. o 

Note that the condition (b) of Statement 2 is automatically satisfied for 
H = 1/2, because y is a stable Levy process. The next theorem shows that this 
condition also holds for general self-similar processes with stationary increments 
(SSSi). It is the case for H = 1/2 again. 

Theorem 3. Let £(t), £(0) = be an SSSi-process for which sample paths 
have only discontinuities of the first kind, a.s. Let M be sup{£(s),s G [0, 1]} 
and G be the leftmost position of M: 

G = mf{t e [0, 1] : sup£(s) = M} 

s->t 
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Then 

1) G has a continuous probability density ip(t) in (0,1), and 



m < ^(ajmax^.-^— |J, Vt,se(0,l), 

i.e., ip = or tp > 0. 

2) if ip ^ 0, then the position Gt of the supremum of £(t) in (—Ta, T(l — a)) 
satisfies the following relation: 



P(|G T |<1) = ^(l + o(l)), T-oo. (18) 

Remark. Theorem 3 shows that the asymptotics of type (18) is due to the 
presence of a nonzero distribution density for the position of the supremum of 
the SSSi process in (0,1). That fact was first pointed out in [15] for the process 
bn(x). The proof of the general case is nearly identical with that given in [15], 
so it is relegated to the Appendix. 

The unilateral IFBM process. 




Figure 2: Distribution of the position of maximum G for the IFBM process. 
Shown are the parts of curves related to the interval (10~ 4 , 10 _1 ) together with 
the corresponding values of H 

The process y(x) — bii(s)ds in (0,1) was generated with the same pa- 
rameters H, T and N as in the bilateral case. The computation is illustrated, 
as above, by the distribution of the position of the maximum G for y(x) with 
H = 0.1 - 0.4 in Fig. 2 (left) and with H = 0.5 - 0.9 in Fig. 2 (right). The 
plots clearly reveal the influence of the atom Po(T) at zero for H > 0.5, which 
impedes to see the exponent = 8(H) in the power law asymptotics of Fq near 
zero. One can nevertheless assert that 0(H) — ► at the cndpoints H = and 1. 
The maximum likelihood estimates of 8 in the interval of G : 10~ 3 — 10~ 2 are 
as follows: 
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H : 0.1 0.2 0.3 0.4 0.5 
6 : .09 .13 .18 .21 .23 
8 : .09 .16 . 21 .24 . 25 

We also list the hypothetical values of 9 : 9o = H(l — H) for comparison 
purposes in the above table. The exact result due to Sinai [18]: 9 = .25 for 
H = 1/2 shows that we can still take the error of 9 equal to 0.03. In that case 
the hypothetical estimates y(l) do not contradict the empirical ones. 

Speaking in terms of the process y(l), which can be obtained by replacing 
the time in y(x) with local time of bn(x), the hypothesis 9 (H) means that the 
position of the maximum Gl in y(l), < I < L has the property 

P{G L <l) = L- H C L , L^oo (19) 

where £^ is a slowly increasing function of L. In the limiting case H — > the 
local time l(x) ~ x, hence y(l) = £1 where £ is a Gaussian variate. Consequently, 
P(Gl < 1) ~ 1/2, H — ► 0, which is consistent with the hypothesis (19). 

It would be natural to expect an analytical dependence of 9 on H for the 
IFBM process. Consequently, the hypothesis 6q{H) can also be extended to 
cover the case H > 1/2. That extrapolation is exact for H = 1, because 
y(x) = £,x 2 /2, where £ is a Gaussian variate, so that P(G < x) = 1/2 for 
H = 1. The last result corresponds to 9(1) = 0. 

The rigorous result guarantees that 9(H) /(l + H) is decreasing for H > 1/2. 

Statement 4. (a) The distribution Fm(x\H) for the maximum of 

j/i (i) = V2H + 2 [ b H (s)ds, < i < 1 (20) 
Jo 

increases with increasing H in the interval (1/2, 1) for any fixed x > 0. 

(b) If Fm(x 1+h \H) = x ^ C(x), ie J. where C(x) is a slowly varying 
function, then 6(H) /(l + H) decreases with increasing H in (1/2, 1). 

Proof. The process (20) differs from y(x) by the normalization i?|yi(l)| 2 = 1. 
Let £ q (x) = yi(x e ) where q = 2H + 2, 6 = q /q, q = 2H + 2 and H > H . 

Since IFBM is self-similar with parameter h = H + 1, one has E\^ q (x)\ 2 — 
\x\*° = E\£ qo (x)\ 2 

We show in the Appendix that 

EZ q (x)Z q (y) > Et qo (x)t qo (y) (21) 
when H > 1/2. In that case the Slepian lemma [13] yields 

P(m&xf„(x) < u) > P(m&xf 0n (x) < u). 
[oi] v [01] ^ yow 
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However, max[oi]£ 9 (x) = maxpi] \J2H + 2 J* bfj(s)ds which proves the first 
part of the statement. The second part is an obvious corollary of the first. 

5. Conclusion 

We were testing the hypothesis that the maximum M of the integral of 
fractional Brownian motion with index H has the distribution Fm(s' 1+h ') = 
x 8 ( H ) £ H (x), x — > in I 3 0, where Lh is a slowly varying function. We have 
presented theoretical arguments and computational evidence to support and 
refine the hypothesis as follows: 9(H) = 1 — H for / = (—1, 1) and 9(H) = 
H(l — H) for I = (0,1). The computational part of the problem faces the 
following difficulties. Due to discrete time (a step of A), the analogue of M in 
the grid case (M) has a nonzero probability P(M = 0) of the same order as 
Fm(A 1+h )- It causes a slow (power law) convergence of the distributions of 
M and M. Further, there exists an interval (0, to(A, H)), t — > as A — > 0, 
where P(M < x 1+H ) has a power law behavior with exponent 0(H) > 9(H), 
H > 1/2. This makes the choice of the interval where 9(H) is to be estimated 
more difficult. When A = 0.00012 and the number of samples is N = 50000, 
the accuracy of 9(H) is ~ 0.03 for the most favorable range of H: H < 1/2. 

6. Appendix 

Proof of Theorem 2. A distribution function (here Fq) is differentiable al- 
most everywhere, because of monotonicity. Suppose this is true for the point 
x . Consider x < xq, A = x /x > 1. The self-similarity of £(x) yields 

P(G(0, 1) G dx) = P(G(0, A) e Xdx) < P(G(0, 1) e \dx) = ip(x )—dx (22) 

x 

Here, G(a,b) is the leftmost position of the supremum of {£(x),x G (a, b)}. 
Consequently, the distribution of G(0, 1) is absolutely continuous in (0,xo). 
Points like xq are dense in (0,1). Consequently, Fc(dx) — ip(x)dx, x G (0, 1). 

In virtue of (22), ip(x)x is a nondecreasing function, i.e., the discontinuities 
in ip arc at most denumerable, while finite limits on the left and the right exist 
at the discontinuity points. The fact that the increments of £(x) are stationary 
yields 

p|g(0,1) e ctej = P{G(-a, 1 - a) e dx - a} < P{G(0, 1 - a) e dx - a} 

for any < a < 1. One has 

, / n , fx — a\ 1 

i)(x)<M- — - — 

\ 1 — a/1 — a 
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at continuous points of ip. Multiply both parts by (1 — x): 

(i - xMx) < V (^) (i - ^) = V(y)(i -y), y = ^<x. 

Combining both inequalities, one gets 

ip(x) < V-(y)max^,i-|^) (23) 

at all points where x and y are continuous. In particular, ip(x + 0) < ip(x — 0) < 
+ 0), i.e., -0 is continuous in (0,1). If ip(xo) — 0, x € (0, 1), then one has 
^(a:) = 0, a; e (0, 1) from (23). Consequently, the following alternative holds: 
cither ip = or ip > in (0,1). The second part of Theorem 2 is an immediate 
corollary of the first part and the self-similarity of £(x), see [15]. 
The proof of (21) in Statement 4- 

Let £q(t) = /„ T b H (s)ds, r = t fl ', = 9o /g < 1, q = 2H + 2. The 
correlation function /3 q (t, s) of £, q (t), can be written as 

2t-*>&(t, s) = -^(/ + P 90 - 8 ) + [(1 - P 6 )* 79 - (1 + P 90 )]-^ (24) 
qo-a qo - V 

where p = s/t. Because f3 q is symmetric in t, s, we put p < 1. We will show 
that /3 q (t,s) > (3 qo (t,s), if q > qo > 3 or, which amounts to the same thing, 
H>H > 1/2. 

We have in virtue of (24): 

2 [/?,(«,«) - (3 qo (t,s)] = («- 9o)(l - P%o - 1) -1 -Ri + <7o(go - l)" 1 ^. 

Here J2i = (1 - y^ 1 - y^q - l)" 1 - (y*- 1 - y^X? - q )-\ y = p e , 
y = 1 — y and 

R 2 = f [p a - p q °- a - (1 - p a ) 9 °-V] da lnl/ p. 

We now are going to show that i?i > 0, and i?2 > 0, if H > 1/2. 
Consider R 2 . Put p Q = u. Since < 6 < a < 1 and < p < 1, it follows 
that < u < 1. The integrand in R 2 becomes 

u[l - u go/a - 2 - (1 - uY ' 1 ] > u[l - u g °- 2 - (1 - u) q °- 2 ] 

> u[l -max(l,2 3 -« )] > 

The last estimate is true, because 3 — qo = 1 — 2H < when H > 1/2. 
Consequently, R 2 > 0. 
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Consider R\. The function R\{y) is positive around and 1: 



f y 2 {q -l)/2 + O{y^), y^O 
Hl -\ y + Oig*- 1 ), y-^0. [2b > 

Consequently, R\ > 0, if the function has a single local extremum in (0,1). 
Let z = (1 - y)- 1 e (1, oo). Then 



^°- 2 ^i?i - [(* - I) 9 " 2 - 1] - [z q - q °(qo ~ 1) - q ~ l](q - go)" 1 := f{z) 
We now show that f(z) has a single root in (1, oo). The function 



f(z] _ f -(* -l)(9o- l) + 0{(z-lf H ), z^l 

■ nz >~ { (z- 1)"- 2 (1 + o(l)), z^oo [2b > 

changes sign in (1, oo). The equation f'(z) = or 



(q-2)(z-iy- 3 = (q Q -l)z*-*>- 1 (27) 

determines the local extremums of f(z) in (1, oo). Two strictly monotone func- 
tions occur in (27): on the left is a function that increases from to oo, because 
q — 3 = 2H — 1 > 0, while the nonnegative function on the right decreases 
toward zero at infinity, because q — q — 1 = 2(H — H ) — 1 < 2 • 1/2 — 1 < 0. In 
that case, however, (27) has the single root 1 < z* < oo. In virtue of (26) z* is 
the point of minimum, where f(z*) < 0. The function f(z) is strictly increasing 
from f(z*) < to oo in (z* , oo), so that the equation f(z) = has a single root, 
as was to be proved. 
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